arxiv 



Quantitative model for anisotropy and reorientation thickness of the magnetic 
moment in thin epitaxial strained metal films 
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A quantitative mathematical model for the critical thickness of strained epitaxial metal films is 
presented, at which the magnetic moment experiences a reorientation from in-plane to perpendicular 
magnetic anisotropy. The model is based on the minimum of the magnetic anisotropy energy with 
respect to the orientation of the magnetic moment of the film. Magnetic anisotropy energies are 
taken as the sum of shape anisotropy, magnetocrystalline anisotropy and magnetoelastic anisotropy, 
the two latter ones being present as constant surface and variable volume contributions. Apart from 
anisotropy materials constants, readily available from literature, only information about the strain 
in the films for the determination of the magnetoelastic anisotropy energy is required. Application 
of continuum elasticity theory allows to express the strain in the film in terms of substrate lattice 
constant and film lattice parameter, and thus to obtain an approximative closed expression for the 
reorientation thickness in terms of lattice mismatch. The model is successfully applied to predict 
the critical thickness with surprising accuracy. 

PACS numbers: 75.70.-i, 75.30.Gw, 75.40.Mg 



I. INTRODUCTION 

The efficiency of magnetic data storage media depends 
critically on the presence of a perpendicular magnetiza- 
tion, this is, the easy magnetic axis is perpendicular to 
the surface of the magnet in the absence of an external 
magnetic field. Generally, the magnetization lies in some 
preferred direction with respect to the crystalline axes or 
to the external shape of the magnet - this is known as 
magnetic anisotropy. The energy involved in rotating the 
magnetization from a direction of low energy (easy axis) 
towards one of high energy (hard axis) is typically of the 
order 10~^ to 10~^ eV/atom and thus a very small correc- 
tion to the total magnetic energy [1] . The physical origin 
for the anisotropy is the symmetry breaking of the rota- 
tional invariance of the dipole-dipole interaction and the 
spin-orbit coupling with respect to the spin quantization 
axis. Atoms at interfaces in ultrathin films or multilay- 
ers show a much stronger anisotropy than atoms in the 
bulk. Also, atoms in strained lattices show a stronger 
anisotropy. 

While it has been an established fact for many years 
that epitaxial strain in ultrathin metal films can have 
significant impact on the stabilization of perpendicular 
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magnetization [2-5], many relevant publications still ig- 
nore magnetoelasticity for the interpretation of magnetic 
anisotropy data. 

The relationship between the magnetic anisotropy and 
elastic strain is also important in structural geology. 
Elastic strain is a marker that permits to draw conclu- 
sions about deformation history of rocks. Since such 
strain data are not always directly available, quantita- 
tive modeling helps to retrieve strain data from magnetic 
anisotropy data [6] . Reviews about the interplay between 
mechanical stress and magnetic anisotropy in thin films 
and in tectonics are given by Sander [7], and by Bor- 
radaile and Henry [8], respectively. 

In the next section, the specific contributions to mag- 
netic anisotropy energy will be given and summarized. 
The magnetoelastic energy will be expressed in terms of 
lattice spacings in film and substrate in order to account 
for the strain in films. The derivative of the total mag- 
netic anisotropy with respect to the film crystal axes, 
i.e. the minimum of this energy will yield a thickness at 
which the magnetization switches from in-plane to out- 
of-plane. The terms out-of-plane and perpendicular are 
used equivocally in this manuscript. The Bain path is a 
functional relationship between in-plane and out-of-plane 
lattice distances of a bulk metal phase and its corre- 
sponding tetragonally distorted equilibrium phases un- 
der epitaxial stress. It can be regarded as a sequence 
of tetragonal states produced by epitaxial stress on an 
equilibrium tetragonal phase [9-11] and serves as a crite- 
rion for the identification of equilibria phases of strained 
structures and for the interpretation of LEED data [12]. 
The concept is being improved still to date and allows, 
for instance, deeper insight in phase transitions during 
epitaxial growth [11], or just the conversion of film thick- 
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ness data from Monolayers to Angstrom [13]. For the 
present model, the Bain path plays a key role. To the 
best of the authors knowledge, the present model is the 
first of its kind which gives a quantitative prediction of 
the critical reorientation thickness. 



II. FORMULATION OF THE MODEL 



A. Magnetocrystalline anisotropy energy 



The magnetocrystalline anisotropy energy GY^ysii^M) 
is written as expansion of the cosine of directions ai, a2, 
and as of Qm- 



with demagnetization field Hd(r) and magnetization 
M(r), depends on the external shape and geometry of 
the magnet. Ultrathin films can be well approximated 
by a plate with infinite lateral extension, such as: 



T^V ■ 2, 

^ shape 



(5) 



Eq. (5) describes the shape anisotropy energy of a sin- 
gle atom in a monolayer of such a film. Considering in- 
finitesimally thin slices yield generally the same angular 
dependence: 
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Qm is the unit vector of magnetization direction with 
components ai, a^, as or the polar angles 9 and (j). For 
crystals with cubic symmetry, the expression in Eq. (1) 
simplifies to 
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with the coordinate system being in line with the crystal 
axes, and: ai=sin^ cos Q:2=sin^sin0, and q:s=cos^. 
The vector of magnetization, M, includes with the sur- 
face normal of the (001) plane the angle 9 and with the 
[010]-direction the angle (j)- The expression in Eq. (2) 
converges well and thus permits to interrupt the expan- 
sion after the third order. Ko, Ki etc. are the magne- 
tocrystalline anisotropy constants. Note, that these con- 
stants are temperature dependent materials parameters. 
But here we will treat the case for ambient temperature 
only. 

Eq. (2) describes the magnetocrystalline anisotropy 
energy or a single atom of one monolayer in the volume 
of the film. This contribution grows proportional with 
the film thickness. Atoms at surfaces and interfaces are 
not described properly by Eq. (2) anymore. Because 
of the reduction of symmetry at surfaces and interfaces, 
second order terms come into play. For a surface with 
tetragonal symmetry, siich as the cubic (001) plane, the 
surface magnetocrystalline of a single atom is 
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C. Magnetoelastic anisotropy energy 

For the strained lattice of a magnetized body, the en- 
ergy terms may depend on the strain tensor e and on Q,m ; 
this is the magneto-elastic energy. For small lattice mis- 
matches, the energy terms can be expanded in spherical 
harmonics and in powers of e: 

Gmagn.eli^M,^) = ^ Bij^k,l^i,jO:iOij + (7) 

i,j,k,l 

The crystalline symmetry manifests itself in a coupling 
of the expansion coefficients. For a cubic system, the 
expression reads 

Gmagn.eli^M , = 

-Bi(eiiaf -I- e22ai + ess^s) 
-I- 2B2(ei2Q!lQ!2 + ^23^2^3 + ^siC^sCii) + (8) 

The €ii describe the strain and the Cij describe the shear 
deformation. 

We now want to derive an expression for the magneto- 
elastic anisotropy energy in terms of lattice deformation. 
Let us consider here a face centered cubic lattice. The 
nearest neighbor (NN) distance of the atoms in a cubic 
lattice with no strain be oq. For the lattice; under strain, 
the NN distance in the tetragonal plane be a. For pseu- 
domorphic growth, the latter value is the NN distance of 
the atoms in the substrate. The expansion is equal to 
the lattice mismatch: 
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For the contraction, we find 



a — ao 
ao 



(9) 



B. Shape anisotropy energy 
The shape anisotropy energy 

Gshapeii^M) = dVm{v) ■ Hd(r) 



(4) 
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The variation of c with a is given by the Bain path, 
this is [9] 
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so that we may write 
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Here, Cn and Ci2 are the elastic constants for the corre- 
sponding film material. Some elastic constants and lat- 
tice parameters for Fe, Co, and Ni are listed in Table 
I. 



Element T(K) 


Cll [GPa] 


ci2 [GPa] 


a [A] 


bcc Fe 520 


230 


135 


2.87 


fee Fe 1428 


154 


122 


3.65 


fee Co 300 


242 


160 


3.54 


hop Co 300 


307 


165 


2.82 


fee Ni 300 


247 


153 


3.524 


bcc Ni 300 






2.78 



Table I: Elastic constants and lattice parameters for some 
Fe, Co, and Ni phases [14, 15], hep Co after [16]. 
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We neglect shear effects in the distorted films and in- 
sert the expressions from Eqs. (9)-(ll) in Eq. (8) and 
obtain in first order approximation 



Si <! ( — - 1 I sin- 
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D. Summary of contributions 

All contributions considered so far contain bulk- and 
surface/interface contributions, except for the magnetoe- 
lastic anisotropy energy, for which no surface- or interface 
contributions were available from literature. For a film 
with d monolayer thickness, the anisotropy energy per 
unit surface is thus 



G = d-{G, 
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With the exception of GY^ysi, wc insert the expressions 
for the particular contributions in Eq. (14) and obtain 
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Omission of GYj.ygf is justified for ultrathin films, because 
the error inferred by this contribution is negligible. With 
increasing thickness, however, this contribution will dom- 
inate together with the shape anisotropy, also, because 
magnetoelastic contributions will become smaller due to 
relief of strain in thikcer films. 

We can replace cos^ by (1 — sin' 0) in Eq. (15), so 
that G{a, ao, 0) becomes a function linear in sin^ 9. It is 
of interest now to know the conditions for which the angle 
6 becomes 0. Then, the magnetic anisotropy energy has 
a minimum, and this state corresponds to perpendicular 
magnetization. To determine, for which particular angles 
6 the magnetic anisotropy energy becomes a minimum, 
we calculate the first derivative of G{a^ ao, 9) in Eq. (15) 
with respect to 9. The derivative vanishes for 0=0° and 
90°. For these cases we obtain an explicit expression for 
the critical thickness, at which we have a reorientation of 
the magnetization: 
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The expression in Eq. (16) thus allows the prediction 
of the critical magnetic reorientation thickness, within 



the limits and approximations made in the model. More 
generally, Eq. (15) allows to determine the direction of 
magnetization in a thin film for any given set of param- 
eters. 



III. RESULTS AND DISCUSSION 

The model is applied to thin films of Fe, Co, and Ni, 
and compared with experimental data taken from liter- 
ature. Structure and magnetism of these three metals 
have been subject of intense research for decades, c.f. 
[7, 17] and a whole body of experimental data is readily 
available in literature. 

Wc begin with Figure 1, which shows a schematic dia- 
gram of the critical reorientation thickness dcrit (Eq. 16) 
as a function of a/oo for four different sets and signs of 
Ks and B. Phase ranges for perpendicular magnetization 
are marked with _L ; ranges of in-plane magnetization are 
marked with ||. The modulus for the anisotropy energies 
was chosen arbitrarily to iff-j^^t = 4.0 x 10""*, B = 0.20 
X 10-'', and K^f^^^^ = 8.0 x lO""' eV/atom. The ratio 
a/tto is a measure for the strain in the film. For a/uo < 
1 , the film in-plane lattice parameter is compressed with 
respect to the substrate. For a/ag > 1, it is expanded. 
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The regions of in-plane and out-of-plane magnetization 
are displayed depending on the lattice mismatch and the 
sign of the anisotropy constants, 
into account here. 



G^ryst taken 



For these systems, the magnetocrystalline anisotropy 
energy is dominant in terms of modulus. The posi- 
tion a/ao=l denotes the relaxed state of the film and 
is marked by a dotted vertical line. In addition, a solid 
vertical line marks the specific ratio a/ao, beyond which 
a polar magnetization is not possible anymore, regardless 
the film thickness. This is because in Eq. (16), the shape 
anisotropy and magnetoelastic anisotropy energy cancel 
out in the denominator. 

P. Bruno [1] has compiled a number of anisotropy con- 
stants, displayed in Table II and III, which are used in 
the present work. 



Ki 

K2 

K3 
K3' 

shape 

Bi 
B2 
B3 
B4 



bcc Fc 
4.02 10-" 
1.44 10"* 
6.60 10-9 



hep Co 



fee Ni 



(a) 
(a) 
(a) 



5.33 10"'^ 
7.31 10-f' 



(b) 
(b) 



-8.63 IQ-" (a) 
3.95 10-6 (a) 
2.38 IQ-^ (a) 
6.90 10-^ (c) 



8.40 10-^ (c) 
-8.86 10-^ (c) -5.85 10"^ (c) -0.74 10"* (c) 



-2.53 10"* 
5.56 10-* 



-5.63 10-* 
-2.02 IQ-^ 
-1.96 10-3 
-2.05 10-3 



6.05 10 
6.97 10 



Table II: Volume contributions of anisotropy energy 
constants in [eV/atom]. Data for magnetocrystalline 
anisotropy energy arc valid for T=4.2 K. (a) [18], (b) 
[19], (c) [20], (e) [1]. Magnetoelastic constants were ob- 
tained from [1] and [21-24] and are valid for 297 K. 

The magnetocrystalline anisotropy energy is of the or- 
der 10-6 eV/atom and is always dominated by the shape 
anisotropy energy by one (Ni) or two (Fe) order of mag- 
nitudes. They have always a negative sign and thus favor 
magnetization in the film plane. If we consider a Ni film 
with a lattice mismatch of f—l%, the negative magne- 
toelastic anisotropy energy of Ni will overcompensate the 
shape anisotropy energy by a factor of ten. The surface 
and interface contributions to the anisotropy energy are 
constants that add up in the overall anisotropy energy 
of the thin film. Table III gives an overview of the sur- 
face/interface anisotropy constants of a number of sys- 
tems studied in the past, mostly one single crystalline 
substrates. 



System 


K« [mJ/m^l 


Reference 


Fe(001)/UHV 


0.96 


[251 

L "J 


Fe(001)/Au 


0.47, 0.40, 0.54 


[251 

L J 


Fe(001)/Cu 


0.62 


[251 

L "J 


Co/Au(lll) 


0.42 


[261 


Co/Cu(lll) 


0.33 


f4. 27] 

L 'J 


Co/Ni(lll) 


0.31, 0.20, 0.22 


[28-30] 


Ni(lll)/UHV 


-0.48 


[31J 


Ni/Au(lll) 


-0.15 


[32] 


Ni(lll)/Cu 


-0.22, -0.3, -0.12 


[32-34] 


Ni/Cu(001) 


-0.23 


[34] 


System 


Kc;,„„p [erg/cm^] 


Reference 


Fe(001)/UHV 


-0.27 


[35] 


Fe/Ag(001) 


-0.12 


[35] 


Ni(001)/UHV 


-0.017 


[36] 


Ni/Cu(001) 


-0.025 


[36] 



Table III: Surface contributions of anisotropy energy con- 
stants [1]. Proper conversion yields the order of approx- 
imately IQ-* eV/atom. 

Ultrathin films often exhibit perpendicular magnetiza- 
tion up to several monolayers, which is stabilized by the 
positive surface and interface anisotropy energy. For Fe 
and Co, elastic strain works against this effect, and there- 
fore perpendicular magnetization is not favored anymore. 
This is not the case for Ni. Here, we always register a 
negative magnetocrystalline surface/interface anisotropy, 
and a perpendicular magnetization is favored. Elastic 
strain, such as in the case of epitaxial growth, causes 
a positive contribution of the magnetoelastic anisotropy 
energy and thus enhances this effect. Since the magne- 
toelastic anisotropy energy is a volume contribution, this 
is, it increases to the same extent as the volume or the 
thickness of the film, thick films of arbitrary thickness 
with perpendicular magnetization can be grown - at least 
theoretically. The practical limitation is that thick films 
will not maintain the necessary lattice mismatch, i.e. the 
elastic strain, and relief most of the elastic energy by 
dislocations [37, 38]. 

In Fig. 2, we show the variation of the critical thickness 
(curved lines) for Fe, Co, and Ni as a function of a/ao. 
Solid lines indicate dcru after Eq. (16), while dashed 
lines take also the volume contribution of the magne- 
tocrystalline anisotropy into account. The vertical solid 
and dashed lines indicate the ranges where the denomi- 
nator in Eq. (16) vanishes and beyond which, according 
to the model, no perpendicular magnetization is possible. 

First, we note that thin Fe films have a perpendicular 
magnetization (_L). Thick Fe films have the magnetiza- 
tion vector in the film plane (||). This is the situation as 
represented by Figure 1, plot c). The critical thickness 
for the reversal of the magnetization vector is increas- 
ing with decreasing strain in the film. According to our 
model, particularly an ultrathin Fc film has a perpendic- 
ular magnetization. The dotted line was obtained after 
Eq. (16). The volume contribution of the magnetocrys- 
talline anisotropy energy was not taken into account. But 
we can estimate the influence of this contribution on the 
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0.8 



0.9 



1 

a/a, 



1.1 



1.2 




0.8 



0.9 







a/an 



1.1 



1.2 



FIG. 1: Schematic diagram of the critical reorientation thickness. Each sub-figure represents a specific case with combinations 
of sign of magnetoelastic anisotropy constant Bi and surface/interface anisotropy constant K^ryst- Shape anisotropy is taken 
negative for aU cases. 



critical film thickness. For a 10 ML thick film, the con- 
tribution does not exceed 4.02 10^'^ eV/atom. The solid 
line in Fig. 2 takes this contribution into account and is 
thus an upper approximation for the critical film thick- 
ness. 

Figure2, plot a), contains four experimentally ob- 
tained values for the critical reorientation thickness of 
bcc-Fe films, i.e. Fe/Cu(001) [39], Fe/Cu3Au(001) [40], 
Fe/Ag(001) [41], and Fe/Pd(001) [12]. All of them are 
located below the curves for the calculated critical thick- 
ness. For the system Fe/Cu(001) and Fe/Cu3Au(001) it 



is believed that the magnetic properties depend on the 
preparation conditions of the films, which could not be 
taken into account in the model. In addition, for the 
latter case, interdiffusion of Au atoms in the Ni film is 
believed to have changed the magnetic properties of the 
Ni film [12]. It also should be mentioned that the tem- 
perature dependence of the anisotropy constants was not 
taken into account in the present model. But it is strik- 
ing that the experimentally obtained critical thickness 
increases with decreasing film strain, as predicted by our 
model, i.e. experimental data and calculations show the 





2: Critical thickness for Fe (plot a), Co (plot b), and Ni (plot c). 
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same trend. One thorough confirmation of the clear sepa- 
ration between surface and volume contributions is given 
by Thomas et al. [42]. They found that the uniaxial 
magnetic anisotropy of very thin Fe/GaAs(001) is not a 
result of magnetoelastic coupling nor a result of shape 
anisotropy, but a result of interface anisotropy. On the 
other hand, they find that the further evolution of mag- 
netic anisotropy with increasing Fe film thickness is a re- 
sult of the competition between magnetoelastic coupling 
and interface anisotropy. It is worthwhile to mention 
that their work was based on a thorough X-ray structure 
analysis of the films. 

In the middle plot in Fig. 2, plot b), the critical thick- 
nesses for Co films with and without correction for the 
volume contribution of the magnetocrystalline anisotropy 
are given. Ultrathin Co films have an easy magnetic axis 
perpendicular to the film plane, and thick planes are mag- 
netized in the film plane. This is a situation similar as 
with Fe films, as represented by Figure 1, plot c). The 
experimental data are taken from Co/Au(lll) [4] and 
Co/Cu(001) [43]. In the former case, elastic strain in 
the film is anticipated, but not quantified. The lattice 
mismatch for this system would be theoretically about 
14%. But it is unlikely that the film would maintain 
such a large mismatch; pseudomorphic growth is thus 
ruled out for the face centered phase of Co. Therefore, 
the data point (3.1 ML) was set to the natural lattice pa- 
rameter of Co (3.544 A). The lattice mismatch of fcc-Co 
on Cu(OOl) is approximately 2%, and the critical thick- 
ness was found to be 3.4 ML. Thus, both data points 
are in the permissible region between the two extreme 
cases of the model. Chen et. al. [44] observe in-plane 
anisotropy for Co/Pt(lll) films larger than 3.5 ML, and 
perpendicular anisotropy for thinner Co films. Capping 
these Co/Pt(lll) films with thick enough Ag layers how- 
ever increased the reorientation thickness to as far as 7 
ML, depending on the thickness of the Ag capping layer. 
Lee et al. [45] found for the systems Co/Pt(lll) and 
Co/Pd(lll) reorientation thickness from perpendicular 
to in-plane of 5 ML and 12 ML, respectively. Sandwiches 
of Pd(lll)/Co/Pd were shown to have perpendicular 
magnetization up to 13 - 15 ML [46], which translates 
to approximately 5 ML. For Co/Cu(lll), a critical reori- 
entation thickness of 5.5 ML is reported [45]. Kohlhepp 
et. al. explain the strong perpendicular anisotropy of 
Co films on Pd(lll) substrates, capped with either Fe or 
Cu, with the alloying of Co and Pd at the Co/Pd(lll) 
interface [47]. 

However, in this thickness range, a gradual phase 
transformation from the fee phase to a hep phase takes 
place [17, 45]. It is also possible that a bee cobalt phase is 
synthesized [48] in a particular thickness range, which is 
not necessarily stable, however [49] . Cobalt was also one 
of the first systems where magnetic domain formation 
was found and related with magnetic anisotropy [50-53] . 

The critical thickness of Ni is shown in the bottom part 
of Fig. 2. Ultrathin Ni films are in-plane magnetized, and 
thick films are out-of-plane magnetized. This situation 



is represented in Fig. 1 by plot b). So far we have not 
yet addressed the differences in the interface anisotropy 
energy, which arises when different chemical elements are 
used as substrates, or even when the surface of a film is 
capped with a different element. 

Table III lists a number of surface and interface 
anisotropy constants for various systems and chemical 
species. Since these anisotropy constants are specifically 
given for the coordination numbers and packing densities 
of the constituents of the interface, including their chem- 
ical species, they cannot simply be adopted for systems 
with different coordination numbers etc. 

The dashed and solid lines in the bottom plot in 
Fig. 2 represent actually four lines for the critical thick- 
ness, each with surface anisotropy constants for either 
Ni/Cu3Au(001) and Ni/Cu(001).' The differences of 
these curves are so minute, however, that the curves can- 
not be distinguished anymore. This holds even for the 
curves which take the magnetocrystalline anisotropy into 
account. The reported critical thickness for Ni/Cu(001) 
is 12 ML [54], in good agreement with our model. 

There are no anisotropy constants data available for 
the system Ni/Cu3Au(001), but we can try to find a rep- 
resentative average for this system, built from the linear 
combination of data for Ni(lll)/UHV, Ni/Au(lll) and 
Ni/Cu(001), which are readily available (Table III). We 
assume that the CusAu substrate surface is occupied to 
3/4 of Cu atoms and 1/4 of Au atoms. For simplicity, we 
ignore here that the CusAu substrate may have a Au-rich 
segregation profile on its surface [55], and that alloying 
may occur between Ni and Cu or Au. 

The Cu(OOl) surface has two atoms per unit mesh and 
a lattice parameter of 3.61 A. The number density on 
the surface is thus 1.535 10^^ atoms/cm^. The inter- 
face anisotropy energy constant for Ni/Cu(001) is thus 
-0.935 10~^ eV/atom. Taking the 3/4 occupation prob- 
ability of Cu on Cu3Au(001) into account, we obtain - 
0.701 10-4 eV/atom. 

Nickel does not grow pseudomorphic on Au(lll) due 
to the large lattice mismatch of 15%. Instead, we as- 
sume that Ni grows on Au(lll) in (001) orientation 
with its natural lattice parameter of 3.5238 A. The 
same procedure yields -0.581 10"'* eV/atom for the 
interface anisotropy energy of Ni/Au(001), or -0.145 
10-"* eV/atom with 1/4 occupation of Au. 

For Ni/Cu3Au(001) we thus obtain an occupation 
weighted value of -0.846 10"^ eV/atom for the interface 
anisotropy energy constant. 

Based on the same considerations, we yield -1.6107 
10~*eV/atom for the surface anisotropy of the Ni film 
vs. the UHV. 

Thus, the surface anisotropy and interface anisotropy 
of Ni/Cu3Au(001) amount together to -2.4572 
10-'*eV/atom. 

Eq. (16) neglects the magnetocrystalline anisotropy 
energy, which has a dependence of the azimuthal angle cj) 
and thus makes mathematical expressions less transpar- 
ent then intended here. 
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The model does not include contributions to the mag- 
netic anisotropy energy that arise from inhomogeneities 
such as magnetic domain walls or steps and kinks on sur- 
faces, or surface roughness. For thicker films, Eq. (16) 
therefore becomes less accurate. However, Eqs. (13)- 
(16) are quite simple and yet elegant expressions which 
make our anisotropy model transparent enough for the- 
oretical and also practical considerations. It is antici- 
pated that these expressions can be implemented with 
reasonable effort in other formulae, or that other formu- 
lae can be implemented in our anisotropy model. For 
instance, the magnctoelastic anisotropy contribution is 
based on the assumption, that the lattice parameters of 
the film material remain unchanged during film growth. 
This assumption does not really hold, because of disloca- 
tion formation in the film. Implementation of the force 
balance between film stress and the tension of disloca- 
tions [37, 38] should improve the present model. Crys- 
tal parameters and anisotropy constants are parameters 
that depend on the temperature. Crystallographic phase 
transitions may take place during temperature changes, 
in particular when ultrathin films are concerned. Dif- 
fusion and interdiffusion processes may be activated at 
slightly elevated temperatures and thus facilitate alloy- 
ing of film atoms with substrate atoms, which has in- 
fluence on the bandstructure of the film. Finally, the 
magnetic anisotropy constants are throughout tempera- 
ture dependent and may even undergo changes of their 
sign, which then manifests in the occurrence of so called 
Hopkinson maxima in magnetic susceptibility measure- 
ments [57]. While it is difficult to account for all these in- 
stances and circumstances, it should be appreciated that 
the Bain path at least partially admits to quantitatively 
model the magnetic anisotropy behaviour of ultrathin or 
very thin films on single crystalline substrates. 

IV. CONCLUSIONS 

A quantitative model for the dependence of the mag- 
netic anisotropy of thin films was established. It is based 



on anisotropy constants, lattice strain, and film thick- 
ness and takes volume and surface contributions into ac- 
count. Information about the lattice strain is directly 
implemented in the model as a function of the lattice 
mismatch via the epitaxial Bain path. At the present 
stage, the model does not take into account contribu- 
tions that arise from magnetic domain walls, steps, or 
surface and interface roughness. Also, effects of tem- 
perature variation (phase transformations, temperature 
dependent change of sign of anisotropy constants, etc.) 
are not accounted for yet. But the open design of the 
present model generally allows implementation of these 
contributions and further refinement, on the cost of con- 
ceptual transparency, however. Minimization of the total 
magnetic anisotropy energy of the film yields a closed and 
simple expression for the critical reorientation thickness 
of the magnetic moment of the film. Experimental data 
of Fe, Co and Ni films are in fair agreement with our 
model. The significance of magnetoelastic effects on the 
magnetic anisotropy of thin films is quantitatively veri- 
fied. To account for the infiuence of magnetoelasticity, 
assessment of elastic strain in the films is indispensable. 
A complete structural characterization of the film is thus 
desirable, which permits to assign the correct crystallo- 
graphic phase as well as information about elastic strain. 
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